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ABSTRACT 

In this paper we present semi-analytic models and Monte-Carlo simulations of QSO 
Lya absorption line systems which originate in gaseous galactic haloes, galaxy discs 
and dark matter (DM) satellites around big central haloes. The aim is to estimate the 
number density per unit redshift of Lya absorption lines related with galaxies and to 
investigate the properties of the predicted galaxy/ absorber systems, such as equivalent 
widths W r , projected distances p, galaxy luminosities Lb, as well as absorber redshifts 
z. It is found that for strong Lya absorption lines (W r > 0.3A) galactic haloes and 
satellites can explain ~ 20 per cent and 40 per cent of the line number density of the 
HST QSO Absorption Line Key Project respectively. The population of DM satellites 
is adopted from numerical simulations by Klypin et al. (1999). If big galaxies indeed 
possess such large numbers of DM satellites and they possess gas, these satellites may 
play an important role for strong Lya lines. However the predicted number density 
of Lyman- limit systems by satellites is ~ 0.1 (per unit redshift), which is four times 
smaller than that by halo clouds. Including galactic haloes, satellites and HI discs 
of spirals, the predicted number density of strong lines can be as much as 60 per 
cent of the HST result. The models can also predict all of the observed Lyman-limit 
systems. For strong lines the average covering factor within 250h~ 1 kpc is estimated to 
be ~ 0.36, which is in good agreement with observations. And the effective absorption 
radius of a galaxy (with unit covering factor) is estimated to be ~ WOh" 1 kpc. There 

exist correlations of W r versus p, Lb and z. The models predict W r oc p~ a L B (\-Yz)~ 1 
with a ~ 0.5, j3 ~ 0.15, 7 ~ 0.5. 

To compare with results of imaging and spectroscopic surveys, we study the se- 
lection effects of selection criteria similar to the surveys. We simulate mock observa- 
tions through known QSO lines-of-sight and find that selection effects can statistically 
tighten the dependence of line width on projected distance. This result confirms pre- 
vious suggestions in the literature. After applying selection criteria, the models can 
predict similar distributions of W r , p, Lb, absolute magnitudes and absorber red- 
shifts to those of imaging and spectroscopic surveys. Finally we find that the total 
redshift interval of present observations (~ 5) is not large enough for the models to 
reveal the real relationships if adopting the selection criteria. An adequate total red- 
shift interval might be ~ 10. This may conciliate contraditious conclusions about the 
anti-correlation of equivalent widths versus projected distances by different authors. 
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1 INTRODUCTION 

The origin and nature of low redshift Lya line absorbers 
are still a matter of debate (cf. Chen et al. 1998, hereafter 
CLWB; Tripp, Lu, & Savage 1998). These absorbers are 
thought to arise either from gaseous galactic haloes/discs 
or from the underdense web-like regions of filaments and 
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sheets in the large scale structure of the cosmic matter. The 
first suggestion seems reasonable because galaxies in general 
seem to possess extended gaseous haloes (Bahcall & Spitzer 
1969) or huge gaseous discs (Maloney 1992; Hoffman et al. 
1993). The second suggestion is drawn from the studies of 
high redshift Lya absorption lines. At high redshift, it is 
widely believed that most Lya absorption lines are trac- 
ers of intergalactic hydrogen as suggested by Sargent et al. 
(1980). A diffuse intergalactic medium (IGM) model for the 
Lya forest was investigated by Bi, Borner, & Chu (1992). 
Several authors investigated the scenario of the Lya forest 
produced by the IGM in the context of cosmological simula- 
tions (e.g., Cen et al. 1994; Petitjean, Miicket, & Kates 1995; 
Miralda-Escude et al. 1996; Hernquist et al. 1996; Haehnelt, 
Steinmetz, & Rauch 1996; Cen & Simcoe 1997; Zhang et al. 
1997, 1998; Theuns et al. 1998; Machacek et al. 1999). It is 
natural to extend this model to low redshift. Some insights 
have been provided by hydrodynamic simulations predict- 
ing absorber properties down to z = 0. For instance, The- 
uns, Leonard, & Efstathiou (1998) find that the observed 
decrease in the rate of evolution of Lya absorption lines 
at z < 2 can be explained by the steep decline in the pho- 
toionizing background resulting from the rapid decline in the 
quasar numbers at low redshift (see also Riediger, Petitjean, 
& Miicket 1998). Dave et al. (1999) find that shocked or ra- 
diatively cooled gas of higher overdensity can give rise to 
the majority of strong Lya lines at every redshift. However 
on the observational side, Lanzetta et al. (1995, hereafter 
LBTW) & CLWB claim that luminous galaxies can account 
for at least about 50 per cent (and even more) of the strong 
Lya absorption lines (W r > 0.3 A) observed by the Hubble 
Space Telescope (HST) Quasar Absorption Line Key Project 
(Bahcall et al. 1993; Bahcall et al. 1996; Weymann et al. 
1998). Although the fraction is still quite uncertain because 
of the unknown galaxy space density (or luminosity func- 
tion), the unknown gas absorption cross section of galaxies 
and the uncertainties of the observed number density of Lya 
absorption systems, there is no doubt in principle that some 
absorbers are physically associated with galaxies especially 
for those with column densities above 10 14 cm -2 . Thus the 
origin of Lya lines at low redshift is still an unsolved prob- 
lem. 

In order to clarify this question, it is particularly im- 
portant to identify galaxies giving rise to QSO absorption 
and to analyse their physical properties. So far, tremendous 
efforts have been made to locate responsible galaxies in QSO 
fields (e.g., Morris et al. 1993; Lanzetta et al. 1995; Bowen et 
al. 1996; Le Brun, Bergeron, & Boisse 1996; van Gorkom et 
al. 1996; CLWB; Bowen, Pettini, & Boyle 1998; Tripp et al. 
1998; Impey, Petry, & Flint 1999). In general, these inves- 
tigations of the physical link between an individual galaxy 
and an individual absorption system aim at answering (1) 
whether absorbers are physically associated with galaxies 
and what the percentage of absorption lines is arising from 
galaxies, (2) whether there is an anti-correlation between 
projected distance (i.e., impact parameter) from the line- 
of-sight (hereafter LOS) to the galaxy centre and the ab- 
sorption rest frame equivalent width (hereafter REW). For 
example, Morris et al. (1993) carried out a redshift survey of 
galaxies in the field of 3C273 and found no galaxies within 
a projected distance of 230/i^q 1 kpc of any of the 12 lines to- 
ward 3C273 with REW exceeding 50 mA. One of their con- 



clusions is that the absorbers are not randomly distributed 
with respect to the galaxies, though the absorber-galaxy cor- 
relation is not as strong as the galaxy-galaxy correlation. 
In a similar program but with different LOS, Stocke et al. 
(1995) and Shull et al. (1996) suggested that most of the 
Lya absorbers are located within large-scale galaxy struc- 
tures. In contrast, using an imaging and spectroscopic sur- 
vey (in the field of HST spectroscopic target QSOs), LBTW 
claimed that the fraction of absorbers arising from galax- 
ies is quite high and there is a distinct anti-correlation of 
REW versus projected distance. CLWB confirmed these re- 
sults with more LOSs and concluded that most galaxies are 
surrounded by extended gaseous envelopes of ~ 170/i _1 kpc 
in radius and that many or most Lya absorption systems 
arise in galaxies. Tripp et al. (1998) reached similar conclu- 
sions, but cautioned that selection effects could artificially 
tighten the anti-correlation, and also that the galaxy survey 
may be incomplete. They pointed out that there could be 
fainter galaxies at smaller projected distance (also suggested 
by Linder 1998; see also Impey et al. 1999) which could 
be revealed in a deeper survey. Moreover, they found some 
missing lines from the CLWB samples and from their LOSs. 
These missing lines would weaken the anti-correlation, if in- 
cluded. Recently, Impey, Petry & Flint (1999) studied Lya 
QSO absorbers in the nearby universe (0 < z < 0.22) based 
on the spectroscopy of ten quasars obtained with the God- 
dard High Resolution Spectrograph (GHRS) of the HST. At 
odds with the results of LBTW & CLWB, they concluded 
that nothing in their data would specifically lead to associate 
absorbers preferentially with haloes of luminous galaxies. 

Another very useful tool for the analysis of physical 
properties of absorbers is provided by observations of the 
intervening absorption in multiple LOSs either from close 
quasar pairs or from multiple, gravitational lensed quasar 
images (see Rauch 1998 for the part of review on this sub- 
ject). For example, double LOSs observations (e.g., Dinshaw 
et al. 1995; Dinshaw et al. 1994; Fang et al. 1996; Petitjean 
et al. 1998) have shown that the absorber size is about hun- 
dreds of kiloparsecs, Rauch et al. (1999) analysed spectra of 
images of a lensed quasar at redshift 3.628 and found some 
low-ionization lines arising from the ISM. 

Given these observational results, it is important to 
build theoretical models to understand the origin of the Lya 
absorbers at low redshift. Some theoretical efforts have been 
made to relate absorption systems with galaxies (e.g., Mo & 
Morris 1994; Mo 1994; Morris & van den Bergh 1994; Mo & 
Miralda-Escude 1996; Linder 1998; Linder 1999). Unfortu- 
nately, even for those Lya absorption lines genuinely arising 
from galaxy haloes, we do not know a priori which part 
of the galaxy gives rise to the absorption. In other words, 
it is unclear whether the absorbing clouds are located in 
the outer regions of the halo as infalling clouds or in the 
rotating disc as interstellar medium clouds or in the satel- 
lites. There are some competing models. Morris and van den 
Bergh (1994) estimated that tidal tails can explain ~ 20 
per cent of the low redshift Lya absorbers, but so far there 
is a lack of detailed models. Mo & Miralda-Escude (1996) 
concluded that gaseous galaxy haloes can account for all ab- 
sorbers with HI column density Nm. > 10 17 cm -2 at redshift 
z < 2. Recently a model in which absorption is due to gas in 
an extended disc was proposed by Linder (1998, 1999), who 
argued that high surface brightness galaxies together with 
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low surface brightness galaxies can account for the majority 
of Lya absorption line systems. However this picture only 
incorporates spiral galaxies while there are some absorbers 
associated with E/SO galaxies (cf. CLWB). This model re- 
quires a large number of low surface brightness galaxies and 
it is unknown if a spiral galaxy can possess a disc extending 
beyond lOO/i -1 kpc. 

In current models of galaxy formation, galaxies are con- 
sidered to possess haloes, discs (for spiral galaxies) and satel- 
lites which can give rise to absorption. Motivated by these 
considerations, we perform Monte-Carlo simulations using 
semi-analytic models with plausible assumptions, given our 
current knowledge about the properties of these compo- 
nents. Our aim is to study: (1) which component is most 
important. (2) whether current observational results can be 
explained by the models. (3) whether selection effects (which 
should be applied when pairing an absorber with a lumi- 
nous galaxy) can tighten the correlations between REW and 
projected distance. (4) what kind of future observations are 
needed to discriminate models and to examine the correla- 
tions. 

The paper is organized as follows. In § [| we describe 
our Monte-Carlo simulation methods and give results. We 
construct simulations with allowed parameters and com- 
pare the predicted line number density (^-) with obser- 
vational results for W r > 0.3A Lya lines, Lyman-limit sys- 
tems, and damped Lya systems. Detailed properties of ab- 
sorbers, correlations of equivalent width versus projected 
distance, galaxy luminosity and redshift, are studied. The 
average covering factor is estimated. In § ^, we study se- 
lection effects. After applying selection criteria, we com- 
pare our results with imaging and spectroscopic surveys. At 
last we make some predictions. A discussion is presented 
in § ^. A summary of the results is given in § ^. Through- 
out the paper, we adopt a dimensionless Hubble constant 
h = //o/(100 kms~ 1 Mpc _1 ). Our presentation is mainly 
based on the ACDM cosmogony (with fio = 0.3, £Ia,o = 
0.7, h = 0.7), results based on the SCDM cosmogony (with 
fio = 1.0, CIa,o = 0.0, h = 0.5) are also discussed for com- 
parison. 



2 THE MONTE-CARLO SIMULATIONS 

We start our Monte-Carlo simulations with galaxy samples 
(whose luminosity distribution is consistent with observa- 
tional luminosity functions) along a QSO LOS. These galax- 
ies are placed along the LOS randomly within a cylinder 
volume in co-moving coordinate (then we assign a redshift 
and a projected distance for each galaxy). For one particular 
galaxy whose halo is characterized by the circular velocity 
(derived from the Tully-Fisher relation or the Faber- Jackson 
law), we model its absorbing components (discs, halo gas 
clouds and satellites) in detail so as to determine its gas 
cross section and cloud properties, such as HI column den- 
sity, temperature and LOS velocity. The total equivalent line 
width is then calculated assuming a Voigt profile for each 
cloud. This procedure produces a catalogue of absorber- 
galaxy pairs with information about the absorbing galaxies 
for further analysis. 



2.1 Galaxy samples 

There will be a number of galaxies intersecting a particular 
LOS with random projected distances to a given QSO at 
redshift z q . We consider those galaxies in a cylinder within 
a radius of R cy = 400/i _1 kpc to the QSO LOS, since m our 
models there is no absorption cloud outside of 400h~ kpc 
and the upper limit of the projected distance in imaging 
surveys of absorbers is less than this radius. The number of 
galaxies at redshift z with interval Az is 

AN g = n c (l + zfnR 2 cy ^Az, < z < z q , (1) 
where 

* = I ( 2 ) 

dz (l + z)H(zY y> 

h(z) = #o[av,o + (i-a\,o-a))(i-M) 2 +a>(i+2) 3 ] 1/2 -(3) 

The co-moving galaxy density n c is obtained by integration 
over the B-band Schechter luminosity function 



4>(Lb)cILi 



t>* (Lb / Lb*) 



-°>„- l b/ l b 



* cLLb/Lb* 



(4) 



where ismin is the minimum B-band luminosity. The lu- 
minosities of these galaxies are selected in such a way that 
their distribution is consistent with the luminosity function 
<KL B ). 

The luminosity functions for different morphological 
types over the range —22 < Mb + 51og/i < —14 are 
as follows (Marzke et al. 1998): (1) For late- type galaxies 
(Spiral), cj>* = 8.0 ± 1.4 x 10" 3 /i 3 Mpc" 3 , a = l.lltg;^, 
M B * = -19.431^8 + 5 log ft, and 0.0067L S ,, < L B < 
10.7Ls». (2) For early-type galaxies (E/S0), <f>* = 4.4±0.8x 
HT 3 h 3 Mpc" 3 , a = 1.00±g;g|, Mb, = -19.37±^? +5 log h, 
and 0.007L_b» < Lb < 11.3Ls*. We do not consider 
Irr/Pec galaxies because they are rare (0* = 0.2 ± 0.08 x 
10 _3 fr 3 Mpc" 3 ) f\ 

The luminosity functions above are valid only at very 
low redshift (z < 0.05). They are derived from the recently 
enlarged Second Southern Sky Redshift Survey (SSRS2). 
Some other determinations give higher normalizations. For 
example, the galaxy luminosity functions from the ESO 
Slice Project (ESP) galaxy redshift survey (Zucca et al. 
1997) is characterized by </>* = 2.0 ± 0.4 x 10" 2 fe 3 Mpc" 3 



-19.611^08 and a 



1.22±° ° 6 7 over the redshift in- 



terval z < 0.3. The ESP luminosity functions are in good 
agreement with those of the AUTOFIB redshift survey (El- 



lis et al. 1996) which are characterized by ■ 



2.45 



+0.37 



10" 2 /i 3 Mpc -3 , a 



-1.16l°;J|, M* 



the redshift interval 0.02 < z < 0.15, 0* = 1.48 
10~ 2 h 3 Mpc' 3 , a 



+0.15 
12 

+0.30 



-19.301^ over 



-1 41+ 012 M - 



-19.65to ll over the 

— 3.55_ 9 nn X 



redshift interval 0.15 < z < 0.35, and 
10- 2 fr 3 Mpc- 3 , a = -1.45±{5;i|, M» = -19.38i^f over 
the redshift interval 0.35 < z < 0.75. The luminosity 
function of Lilly et al. (1995) is characterized by 0* = 



* Only normal galaxies are considered in the models. In reality 
however, galaxies may contain HI tidal tails and could give rise 
to absorption. We will discuss the problem in § W 
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2.72 ± 0.4 x 10~ 2 /i 3 Mpc -3 , which is about twice that of 
SSRS2 (Marzke et al. 1998), and a = 1.03 over the redshift 
interval 0.2 < z < 0.5. 

Since the faint-end slope a and the characteristic mag- 
nitude M» of various luminosity functions are not much 
different, we apply the AUTOFIB luminosity function nor- 
malization over 0.02 < z < 0.15 in our simulations and 
assign it to the two morphological types with the same 
ratio for spirals and E/S0 galaxies as in the SSRS2 lu- 
minosity function. Namely, the characteristic luminosity is 
<f>* ~ 1.58 x 10~ 2 /i 3 Mpc~ 3 for spiral galaxies and 4>* ~ 
0.87 x 10~ 2 /i 3 Mpc~ 3 for E/S0 galaxies. We will discuss re- 
sults for other luminosity functions. 

As the luminosity of a galaxy in the cone along the 
LOS and its redshift z g are known, we calculate the appar- 
ent magnitude of the galaxy applying the k-correction and 
cosmological dimming, 



m B = M B + 51og(D L ) + 25 + k B (z s ), 



(5) 



where Dl = (1 + z)Dm is luminosity distance of the galaxy 
in Mpc, and Dm, the proper motion distance is 



Dm = c 



H(zY 1 dz 



(6) 



for a flat universe (Qk = 1 — flo — = 0), which is the case 
in this paper. We adopt B-band k corrections for galaxies of 
different morphological types as in Pence (1976). 

At any given epoch, haloes can be parameterized by 
their circular velocity V c i r , which is simply related with 
galaxy morphological type and luminosity. Empirical rela- 
tions are known between B-band magnitude and LOS veloc- 
ity dispersion, a, of the matter in galaxies both for ellipticals 
and for spirals (Faber and Jackson 1976; Tully and Fisher 
1977). 

The Faber-Jackson relation is 
- M B + 5 log h = (19.39 ± 0.07) + 10(log a - 2.3) (7) 
for ellipticals, and 

-Mb + 5 log ft = (19.75 ± 0.07) + 10(log a - 2.3) (8) 

for SO's (Fukugita & Turner 1991). The circular velocity 
of the halo for elliptical and SO galaxies is V c ir = y/2a. 
Using the E/S0 type luminosity function of Marzke et al. 
(1998), we get 81.6 km s" 1 < V cir < 514.6 km s" 1 (V c * r ~ 
280.1 km s~ ) for Elliptical galaxies and 75.1 km s _1 < 
V^r < 473.6 km s" 1 (V c * r ~ 258.5 km s" 1 ) for SO galaxies. 

For spirals the Tully-Fisher relation is used to derive 
the LOS velocity width Av. We take (Fukugita & Turner 
1991) 

-M B +51og/i= (19.18±0.10)+(6.56±0.48)(log Aw-2.5).(9) 

The halo circular velocity for a spiral is V c ii = An/2. We get 
25.7 km s" 1 < Vox < 425.4 km s^ 1 (V*, ~ 172.6 km s" 1 ) 
using the spiral-type luminosity function of Marzke et al. 
(1998). Note that the upper limits of the circular velocity 
are very large, there are, however, no such large galaxies in 
the sample due to the exponential cutoff in the Schechter 
luminosity function. 



2.2 Gaseous galactic haloes 

We model the gaseous galactic haloes following the work by 
Mo & Miralda-Escude (1996). In such semi-analytic mod- 
els, it is assumed that the gas in a halo has a two-phase (a 
hot phase and a cold phase) structure which, in principle is 
described by the density profiles and the temperature pro- 
file. The density profiles are characterized by the so-called 
cooling radius and virial radius. And the temperature pro- 
file of the hot gas is characterized by the virial temperature. 
Our modeling is summarized as follows (see Mo & Miralda- 
Escude 1996 for more details): 

In cooling flow models, when the gravitational potential 
is important, the core radius of the hot gas profile is similar 
to the cooling radius (Waxman & Miralda-Escude 1995). 
Thus, a self-similar density profile for hot gas is assumed as, 



Ph(r) = p h {r c ) 



where 



2ri 



r(r + r c 



Ph(r c ) 



5/j,kT„ 
2A(T v )t M ' 
We assume 

fgVl 



ph{r c ) = 



47rGr 2 



(10) 



(11) 



(12) 



so that the density of the hot gas at this radius is a fraction 
f g of the total density of the halo. A(T V ) is the cooling rate 
of the gas at the virial temperature. The gas mass fraction f g 
is assumed ~ 0.03 — 0.05. p is the average mass per particle, 
which is ~ 0.6mH with mu being the mass of hydrogen 
nucleus. The cooling radius r c is determined by eq. (11) and 
eq. (12), 



f g A(T v )t M 
5nGp 2 



124.6 kpc V / A-23iM/10Gyrs, 



(13) 



where A_23 is the cooling rate in units of 10 _23 erg s _1 cm 3 . 
The hot gas is taken to be isothermal, so that Th(r) = T V — 
pV^ ir /2k. t M = £/(l + flo), is the time interval between ma- 
jor mergers, since the gas is then shock heated to a stage 
from which it starts cooling. The cooling function A(T„) is 
adopted from Sutherland & Dopita (1993). The age of a halo 
at redshift zi is an integration of eq. (2) 



t = 



(l + z)H(z) 



dz. 



The virial radius is 
r v = V cit /[WH{z)]. 



(14) 



(15) 



When the hot gas is shock heated and starts to cool, it 
will sink to the galaxy centre with velocity u = — fu(r). The 
cooling flow can be described as, 



dp c 
dt 



+ V • (p c u) 



A(T h ) 



■Ph(r). 



We assume u(r) to be a constant: 
u(r) — v c 



(16) 



(17) 



where v c is the infall velocity which must be of the order 
of Kir- In eq. (13), r c is a function of t (here t stands for 
t M , and we omit the subscript hereafter), and we have r c = 
^jjr = r c /2t if we assume a constant T v . Let us set 
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Equation (16) has only one variable r and can be simplified 
to 



da; x 1 dx x 2 (l + x) 2 



(18) 



This equation can be solved analytically (see Appendix A 
for details). 

When r c > r v , r c is simply a parameter rather than 
a physical cooling radius. Thus we assume that the resid- 
ual hot gas at r v is still at the virial temperature and has 
a density such that its cooling time is equal to t, so that 
T h (r v ) = T v and p h (r v ) = (5fikT v )/[2A(T v )t]. In such a 
case, we replace r c in eq. (10), (11) by r v . And we can not 
use eq. (16) to describe the cooling flow. In this case, most 
of the accreted gas will have cooled. Since the total gas mass 
accreted in a halo with circular velocity Kir is fgV 2 ir r v /G, 
the total mass of gas that has been in the cold phase can be 
written as 



M 



/flKir 7 ''-' 

G 



4nx ph(x)dx, 



(19) 



where ph(x) is the density of gas in the hot phase as dis- 
cussed above. The cold gas should form clouds that will 
fall through the halo. Assume a mass flow rate M = M/t, 
and assume that the clouds move to the halo centre with a 
constant velocity v c (which is of the same order as the halo 
circular velocity). Assuming also spherical symmetry for the 
gas distribution, we can write the density of the cold gas as 
(see Mo & Miralda-Escude 1996 for more details), 



Pc(r) = 



f g vl 



4nGr 2 v c t 



mm 
~^2~ 



x 2 ph{x)dx 



(20) 



where r min = min(r„, r c ). 

The cold gas is assumed to be, for simplicity, in spher- 
ical clouds with masses constrainted by various physical 
processes, such as gravitational instability, evaporation by 
hot gas, hydrodynamic instability, etc. Too large clouds will 
eventually collapse to form stars and small clouds will evap- 
orate by heat conduction and also be disrupted by hydro- 
dynamic instabilities. The net effects of these processes is 
to preferentially destroy low mass clouds, so it is possible to 
end-up with a log-normal mass function, like the mass func- 
tion of star clusters observed in the galaxy (e.g., Gnedin & 
Ostriker 1997) , if we begin with a power-law mass function. 
For this reason and for lack of knowledge of the mass distri- 
bution of cold clouds, we assume a log-normal distribution 
of cloud masses M c , 



p(M c )dM c = 



1 



-exp 



ln 2 (M c /M) 

2^ 



dM c 

~m7 



(21) 



Here M is the cloud mean mass. As discussed by Mo & 
Miralda-Escude (1996), the mean mass of the clouds is ap- 
proximately 10 5 — 10 6 Mq . In this paper, we choose M = sev- 
eral xlO 5 M0 and a M = 0.1 ~ 0.3. We also use a constant 
cloud mass in simulations, but the results do not change 
much. 

We model the clouds as spheres of uniform, isother- 
mal photoionized gas confined by the pressure of the hot 



medium, so that pdoudTe = phTh, where T c , the temper- 
ature of the clouds, is about 2 x 10 4 A'. The cloud radius 
is R c = (3A/ c /47rp c i oud ) 1/3 (typically R c ~ 1 - lOkpc 
at a radius within 100 kpc, cf. Mo 1994). The total hy- 
drogen column density through the cloud centre will be 
iVo(H) = i? c p c i ou d/2.3/i, and the H number density is 
n(H) = p c ioud/2.3/x. 

We assume that the cloud is almost completely pho- 
toionized by a constant UV background, and in ionization 
equilibrium. The fraction of hydrogen in the neutral state 
(HI atom), is determined by the flux of the UV background 
ionization field J(v) and nn- We take 

J(v) = J_ 2 i(z)xl0~ 21 (—) a eHergcm- 2 sr" 1 Hz- 1 s" 1 ,(22) 

where v a i is the hydrogen Lyman limit frequency, J-21 (z) = 
0.5 for 2 > 2, and J_ 2 i = 0.5 x [(l + ,z)/3] 2 for z < 2. A break 
in the spectrum at 1/4 = 4Ry (due to continuum absorption 
by He II), with Q(u < V4) = 1 and ©4 = &{v > v A ) = 0.1, is 
included (cf. Miralda-Escude & Ostriker 1990; Madau 1992). 
We take a = 0.5. 

For A(H) < 10 19 cm~ 2 , which is the case for most 
clouds, the cloud is optically thin to the ionizing field with 
ionization parameter U = ^^7, where the ionizing photon 

flux <3?(H) = df, h is Planck's constant. The neu- 

tral hydrogen column density iVm can be derived from the 
code CLOUDY 90 (Ferland 1996). Then we obtain the HI 
column density of a cloud at a distance to the LOS I, 



Nm = Nm{0)y/l-P/R 2 c 



(23) 



where Ahi(0) is the HI column density through the cloud 
centre. 

There might be one or more absorbing clouds in the 
LOS with different velocities with respect to the galaxy cen- 
tre. The velocity structure follows eq. (17) and the LOS 
velocity can be calculated easily. 



2.3 Dark matter satellites 

According to the hierarchical clustering scenario, galaxies 
are assembled by merging and accretion of numerous dark 
matter satellites of different sizes and masses. As pointed 
out by Klypin et al. (1999), this ongoing process does not 
destroy all the accreted satellites. Their paper gives re- 
sults of satellite population around a big galaxy-size halo 
by high-resolution cosmological simulations. The VDF (ve- 
locity distribution function) of satellites within 200/i~ kpc 
and 400ft" 1 kpc is 

n(> U cir , sat ) » 5000 (Kdx.B.t/10 km s" 1 )" 2 75 ^ Mpc" 3 (24) 
and 

n(> Kir.sat) w 1200 (Kir.sat/lO km s" 1 )" 2 75 /! 3 Mpc" 3 , (25) 

respectively, where Kir, sat = (10 — 70) km s" 1 . This number 
of satellites is roughly proportional to (Kir/220 km s -1 ) 3 . 
The velocity dispersion of the satellites is of the order of the 
circular velocity of the central halo. The number of satellites 
in the models and in the Local Group agrees well for massive 
satellites with Kir > 50 km s _1 , but disagrees by a factor 
of ten for low mass satellites with Kir about 10 km s _1 — 
30 km s" 1 (see Klypin et al. 1999 for discusion). Possibly, 
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most of these low mass satellites are dark matter mini-haloes 
or analogy of high-velocity clouds (HVCs) at distance > 
100 kpc in the halo of Milky Way galaxy. In addition, the 
distant HVCs are interpreted as gas contained within DM 
'minihalos'. (e.g., Blitz et al. 1999). It is possible that these 
satellites possess gas but have little or no star formation, so 
that they are faint and can not be found in optical surveys. 
The simulation of the survival of DM satellites has included 
dynamics friction and tidal stripping. Gas in DM satellite 
also suffers from ram-pressure stripping by hot gas in the 
central halo. If the surviving DM satellites can accrete gas 
and the gas is not stripped away, they can contribute to 
QSO absorption because there could be some fraction of gas 
in the neutral state with detectable HI column density (Mo 
& Morris 1994), 



N m « 6 x 10 14 cm" 2 



(Vcir,sat \ 
30 kms- 1 / 



m s ( 



0.05 



-3/4 



J- 



(26) 



where R is the galaxy projected distance (distance from 
satellite centre to LOS) and T4.5 = T/W 45 K. We choose 
T4.5 = 1.0 here. The spatial distribution of these satellites 
in the vicinity of the central galaxy is assumed to follow an 
inverse square law of distance to the centre. 

If the population of DM satellites around big central 
haloes is not predicted correctly by N-body simulations or 
these satellites do not possess much gas (for example, due 
to tidal striping, ram-pressure striping, photoevaporation, 
or supernova-driven ejection), the absorption by satellites 
will be overestimated. 



2.4 Galaxy discs 

Observations of low redshift damped Lya systems show that 
some of these systems are possibly not in normal disc galax- 
ies and their host galaxies are ambiguous (e.g., they could 
be low surface brightness galaxies or faint dwarf galaxies, or 
failed galaxies which are not detected, see Steidel et al. 1994, 
Le Brun et al. 1997, Rao & Turnshek 1998). But in general, 
some damped Lya systems must arise from galaxy discs (e.g, 
Prochaska & Wolfe 1997, 1998 and references therein). In 
this paper we only simulate damped systems arising from 
galaxy discs. 

Unfortunately the HI extent of discs is uncertain so far. 
Several studies involving 21 cm mapping of galaxy discs have 
found 'sharp edges' where the HI column density falls off 
dramatically from a few times 10 19 cm -2 to an undetectable 
level (~ 4 x 10 18 cm -2 ). Such edges have been explained 
by models where the ionizing level increases rapidly from 
the inner optically thick to the outer optically thin regime 
(Maloney 1993; Corbelli & Salpeter 1993; Dove & Shull 
1994a). Maloney (1993) assumed an exponential hydrogen 
disc and used a transition region model to calculate the HI 
column density of NGC 3198. His results are in good agree- 
ment with observations. Other authors have suggested an 
extended power-law disc in the highly ionized regime (Hoff- 
man et al. 1993; Linder 1998; Linder 1999). However we 
will take the plausible assumption of an exponential disc ex- 
tending from the centre to the outer part of a galaxy. For an 
exponential disc, we adopt the model of Mo, Mao, & White 



(1998, hereafter MMW). The galaxy disc is assumed to be 
thin, to be in centrifugal balance, and to have an exponential 
surface density profile, 



E(7?) = Y, exp{-R/R d 



(27) 



Here Rd, £0 and R are the disc scalelength, central surface 
density and distance to the centre respectively. Following 
MMW, we have 



Rd 



and 



8 - 8/l_lkpc (c4)( 



250 km s- 



H(z) 
H 



(— ) 



(28) 



Eo^xlO^cm- 2 ^) (-*-)- 



( ) 

V250 kms" 1 / 



H(z) 



H 



3d 



(29) 



where rrid and jd are the fixed ratios of disc mass to halo total 
mass and disk angular momentum of halo total angular mo- 
mentum respectively. Generally, we choose m d ~ jd ~ 0.05 
throughout this work without considering the instability of 
galaxy discs and evolution of these two parameters. A is de- 
fined as the halo spin parameter, whose distribution is 



p(\)d\ 



2tto x 



exp 



In 2 (A/A) 



<2o\ 



dX 
A ' 



(30) 



where A = 0.05 and a x = 0.5 (MMW). 

A galaxy disc is thought to have a vertical structure. 
It is a good assumption (expect for very flattened haloes) 
to ignore the change in the halo density with a height Z 
above the midplane (cf. Maloney 1993). Then in the limit of 
negligible self-gravity the vertical profile of the gas will be a 
Gaussian, 



n H (R,Z) = nn(R,0)e~ 



l /0~ 2 



where the scale height is given by 
°h(R) 



P CTzz 
Va' 



(31) 



(32) 



Here we assume the core radius of the halo to be much 
smaller than R (cf. Maloney 1992) . The asymptotic velocity 
Va is assumed to be of the same order as the halo circu- 
lar velocity V c i r - And we take the typical velocity dispersion 
a zz ~ 6 km s _ . The midplane density is 



nn(R,0) 



(2tt)VW 



(33) 



where the total hydrogen column density N^ t (R) = 
E(E)/mH. The incident ionizing photons come from the top 



of the gas disc with a flux 



photons cm s . This 



photon flux will ionize the gas to a depth Zi at which the 
column recombination rate equals the ionizing photon flux, 



i.e. 



n 2 l (R,0)e 



z 2 / a 2 1 

hdZ — -<fii,ex- 



(34) 



Here 



Qrcc = 4.I8 x 10 13 cm 3 s" 1 T-°- 72 is the re- 
combination coefficient at a temperature of T e ,4 = 
r e /10,000iY. We assume T e = 20, 000K. In the op- 
tically thick regime, the UV ionizing field is incident 
from one side so that the ionizing photon flux is j^i.ex, 
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while <fii t ex = 5.4 x 10 4 /i y photons cm 2 s 1 and I\ y 
J_2i/(0.04erg cm" 2 sr" 1 Hz" 1 s" 1 ). We define 



fjTa r ecr%(R,0)a h 



1.36 x 10 



3 / Vh 

kpc 



Te®' 72 !^ 1 /Nig, 



where N 18 = A^°7(10 18 cm" 2 ). Then the equation (g4J) be- 
comes 



1 - ¥xi{Zi/<j h ) 



(35) 



where Erf is the error function. 

When b = 1, Zi — 0, one can get a critical column 
density below which the hydrogen will be highly ionized. 
This critical column density is ~ a few x 10 19 cm" 2 (see 
Maloney 1993). 

When b < 1, the disc is optically thick and one can 
derive a Zi from equation ( |35| ) . The HI in disc has a sandwich 
structure. If Z < Zi, all the hydrogen is assumed to be 
neutral (the central HI layer) so that the HI fraction is 



Xhi 



(36) 



Above the central HI layer, which is the case for Z > Zi, 
the hydrogen is highly ionized. Thus assuming ionization 
equilibrium, the fraction of HI is 



XHi^2n H (a + 3)T- ' 72 /- 1 , 



(37) 



given the ionizing field is incident from one side of the disc. 

When b > 1, the whole disc becomes optically thin and 
thus the fraction of HI can be estimated as (given the ion- 
izing field is incident from both sides of the disk) 

Xm -tih^ + S)?;- ' 72 /- 1 (38) 

assuming ionization equilibrium (cf. Maloney 1992). Here a 
is the spectral index of the UV background ionizing field, 
which is ~ 0.5 — 1.5. 

Once Xm is known, we can calculate the total HI column 
density along a LOS by 



N m = / Xmn B {R,Z')dZ', 

J LOS 



(39) 



where the LOS has mid-plane distance Ro, inclination angle 
7 and orientation angle a. For a point (R, Z') along the 
LOS, we have 



R(Z') = y/Rl + (Z'tan7) 2 + 2R Z' tan 7 cos a. 



(40) 



Thus we can calculate iVm numerically. We also do calcula- 
tions using Cloudy 90 (Ferland 1996) and find the result is 
consistent with our calculation of the above sandwich struc- 
ture of the disc using the analytic methods. 

In the optically thick regime, our predicted HI column 
density is in good agreement with the calculation by Mal- 
oney (1993) and the fraction of hydrogen in the HI phase 
is about 2/3. However as we shall see in the next section, 
the gas cross section of this regime could be too large and 
the models predict too many damped Lya systems (DLAs) 
at low redshift. Thus we assume the fraction of the hydro- 
gen column density in the HI phase to be ft in this regime 
and simulate for different values of k to compare with the 
observational number density of DLAs. 



2.5 Rest frame equivalent width 

We 'observe' haloes, discs, and satellites by random LOS 
and obtain HI clouds which contribute to Lya absorption. 
Each absorbing cloud (halo cloud, disc, satellite) is assumed 
to have the Voigt profile. The optical depth of a line is 



2.65 x W-tfrkNj^Vjk) 



(41) 



for hi/jk <C kT, where fjk is the oscillator strength of the 
line and Nj is column density, (f> is the Voigt profile, j, k 
are the lower and upper energy level indexs respectively (for 
Lya line, j — l,k = 2). The REW, defined in frequency 
units is, 



W= (1 



)dv. 



(42) 



In accordance with observational usage, W is defined in 
wavelength units, so the value must be multiplied by \/u. 
For the HI Lya line, A = 1215. 670A, and f jk = 0.4162. 

Because there may be n components in a LOS (a 
direct sum or a blend of two or more lines), we com- 

rt 

pute t(v) — S T Ul i and then calculate W numerically. For 

i=l 

Nhi > 10 19 cm~ 2 , which is the case when LOS intersects 
an optically thick disc, the REW is determined from the 
column density accurately as (Petitjean 1998) 



W r = yjN m /(1.88 x 10 18 cm- 2 )A. 



2.6 The predicted line number densities 



(43) 



There are two kinds of basic observational facts for Lya ab- 
sorbers at low redshift. From spectroscopic surveys in QSO 
spectra one can derive the line number densities per unit 
redshift interval, (^) (for absorbers with W r > 0.3A, for 
metal absorption systems, for LLSs, DLAs). Another kind 
of observation is the imaging and spectroscopic survey in 
the QSO fields from which one can relate absorbers with 
luminous galaxies. In this part of the paper, we compare 
the predicted (4^) with that observed to set constraints on 
model parameters and absorbing components. The only se- 
lection criterion for ( ^f - ) is the lower limit of the line width 
(or HI column densities). However, to compare our results 
with results of imaging and spectroscopic surveys, it is nec- 
essary to study selection effects (such as limitations on ap- 
parent magnitude, velocity separation, angular separation) 
when relating absorbers with luminous galaxies and study- 
ing the properties of absorber-galaxy pairs. This study will 
be presented in § ^. 

Several models are simulated with different absorbing 
components and different parameters in plausible ranges. 
The parameters used to describe the absorbing components 
are as follows: f B (the baryon fraction), Z (the metallicity, 
in unit of Zq), M$ (the mean mass of halo clouds, in units 
of 1O 5 M0), v ln _f (the infall velocity of halo clouds, scaled by 
the halo circular velocity, Kir), ft (the HI fraction of total 
hydrogen in the optically thick part of the galaxy disc), the 
velocity distribution function (VDF) chosen for satellites [we 
denote the VDF of satellites following eq.(^) as VDF1 and 
VDF following eq.(||) as VDF2], and the flux of the UV 
background ionizing field J- 2 i(z) = J° 21 (l + z) 2 (J° 

21 ls 

the flux at 2 = 0). 
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Table 1. Results of disc-absorbers with different re. See text for 
discusion. 



0.1 
0.2 

0.3 
0.4 



,dN\ 
(-)dla 

0.03 

0.06 

0.07 

0.09 



0.17 
0.13 
0.12 
0.09 



> 0.3A) m)u 



0.27 
0.29 
0.29 
0.26 



dz 

0.59 
0.60 
0.57 
0.59 



*(^j-)disc is the line number density for all disk absorbers with 
AThi down to 10 12 cm -2 . 

Five types of models are considered: (1) Model A: discs 
only, (2) Model B: halo clouds only, (3) Model C: satellites 
only, (4) Model D: disks and halo clouds, (5) Model F: disks, 
halo clouds and satellites. Model A, B, C are used to pre- 
dict the fraction of absorption by galaxy discs, halo clouds 
and satellites respectively. Model D and F are used to see 
what fraction the plausible combination of absorbing compo- 
nents can explain observational line number densities. The 
standard model is Model F3 which includes all absorption 
components and has standard parameters: 
f g = 0.05, Z = O.3Z ,M 5 = 5,v inl = 1.0, J° 21 = 0.056, k = 
0.1, and VDF2 (the reasons of choosing these values are 
given below). 

Model A has only one parameter, re (see below for its 
meaning) which is listed in Table [[]. The detailed parame- 
ters of the other models are listed in the second column of 
Table ^. In the table, we only list those non-standard param- 
eters and assume a ACDM cosmogony. But we also discuss 
the SCDM cosmogony in some cases (see notation in the 
table) . 

To get statistically stable results, we 'observe' through 
sufficient numbers of LOSs (typically several hundreds). The 
redshift interval for each LOSs is from to 1, and contains 
approximately 300 galaxies in the simulations with V c ir > 
30kms _1 within a column with a radius of 400 hT 1 kpc 
along a LOS. We simulate 500 LOSs for mode A and 200 
LOSs for model B to F (i.e., with ~ 6 x 10 4 galaxies in total) 
to get stable results. The results are listed in Table |l] and 
Table |^. The details of the models are as follows: 

(i) Model A: As mentioned above, the HI fraction of the 
disc model in the inner part of the disk is ~ 2/3. This pre- 
dicts too many low redshift DLAs (see below and Table |l|). 
However some gas in the disc may form stars or may be in 
molecular phase, only a fraction of the gas in the inner part 
of disk is in HI phase. So we define this fraction as a free 
parameter re, and simulate for a set of re to get constraints on 
it by comparing the predicted DLA number densities with 
the observational results. Note that we adopt the traditional 
definition of DLA with N m > 2 x 10 20 cm" 2 . The observa- 
tional number density of DLAs derived by Rao et al. (1995) 
is n OLA (z) = dN OLA /dz = (0.015 ± 0.004)(1 + 2 ) 2 - 27 ±0- 2 5. 
For a mean redshift of z = 0.5, n DLA (0.5) ~ 0.038 ± 0.014. 
The result of the HST QSO Absorption Line Key Project 
(Jannuzi et al. 1998) at z = 0.58 is about 0.020. Our sim- 
ulation results for 500 LOSs with total redshift interval of 
500 are listed in Table |l| Our predicted number densities of 
DLAs are ~ 0.04 — 0.06, comparable to the observed coun- 
terparts, if re = 0.1 — 0.2. Models with larger re will predicted 
too many low redshift DLAs. Thus a reasonable value of re 
is 0.1 and we use this as the standard parameter. 



(ii) Model B: Using model B, we perform simulations for 
different values of f g , metallicity Z, mean cloud mass as 
well as cloud infall velocity. We only take into account cold 
gas inside r m i n of haloes since we assume that there is no 
cold gas outside this radius. The changing pattern of total 
absorber number with different parameters is conceivable. 
For comparison (see Table ^), the results of models Bl, B2 
and B3 show, that the predicted number density will de- 
cline for larger mean cloud masses because there are fewer 
clouds and thus the covering factor is smaller. The mean 
cloud mass will be several x10 5 Mq. Cold clouds of such 
mass can survive from evaporation by heat conduction and 
gravitational instability, as well as hydrodynamic instabil- 
ity (Mo & Miralda-Escude 1996). The results of models B3 
and B4 show that the predicted number densities do not 
change significantly for Z = 0.3Zq and Z = O.IZq. We use 
Z ~ 0.3Zq hereafter. The results of models B2 and B5 show, 
that the predicted number density will increase when the gas 
fraction f g increases (thus increase the total gas cross sec- 
tion). Comparing B2 and B6, if the infall velocities of the 
clouds is only some fraction of the halo circular velocity Kir, 
the predicted number density increases because there should 
be more cold ga s nowadays in the haloes, which is appar- 
ent from eq. (ECf) . The results of model B7 shows, if the UV 
background ionizing field were lower at z < 1 due to the 
declining quasar density, more absorbers can be expected. 

(iii) Model C: Using model C, we compare results for dif- 
ferent f g and VDF of satellites. The simulation results are as 
follows (see Table |^): model CI (C3) predicts more absorbers 
than model C2 (C4) because of the larger gas fraction [see 
eq. (^)]; model C3 (C4) can produce more absorbers than 
model CI (C2) because more satellites are included. Model 
C3 can predict ~ 40 per cent of the line number density 
of the observed strong Lya lines. Our population of satel- 
lites is chosen from Klypin et al. (1999). If it is true that 
a big galaxy possesses a large number of satellites and they 
possess gas, theses satellites may play an important role for 
strong Lya absorption lines (with W r > 0.3A). 

(iv) Model D: The predicted number density of W r > 
0.3A lines in model D is similar to those of model B since 
the contribution by galaxy discs is only a small fraction (see 
Table ^ & Table ^ . Galaxy halo clouds together with galaxy 
discs cannot explain the majority of observational number 
densities for the strong lines. 

(v) Model F: All possible components are considered in 
this type of model. The predicted absorber number densities 
are larger than those in halo-only models and satellite-only 
models (see Table ^). For example, model F3 predicts that 
the number density of absorbers can be as high as 55 ± 22 
per cent of the observed number density of the HST QSO 
Absorption Line Key Project. Model F5 can predict 59 ± 23 
per cent of the observed strong absorption lines. This is in 
good agreement with the results of LBTW and CLWB. It 
is easy to understand that the predicted number density 
should be higher if the UV background is weaker: For exam- 
ple, in model F6, whose UV background is assumed to be five 
times lower, the predicted fraction can reach 92 per cent, be- 
cause more gas in the galactic haloes is in the HI phase and 
the total absorption cross section of satellite-halos becomes 
larger. This model, however, predicts too many Lyman-limit 
systems (see below). 

(vi) If we choose a SCDM cosmogony, for example in 
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model B8, D4, F7, the total absorber number density will de- 
crease for two reasons. One reason is, that for a fixed baryon 
fraction, the cosmic time becomes shorter than in a ACDM 
model so that less gas can cool (the cooling radius decreases, 
as does absorbing gas cross section). Another reason is that 
in a SCDM universe there are fewer galaxies (about 221) 
along a LOS than in a ACDM universe. In summary adopt- 
ing a SCDM cosmogony reduces the fraction of absorber 
number density to only about 28 per cent in model F7, which 
is too small compared with the HST result. We conclude that 
models with a SCDM cosmogony can not predict sufficient 
absorbers at low redshift under the same assumption. 

The observed (^p) of the HST QSO Absorption Line 
Key Project (Bahcall et al. 1996) for strong Lya lines (W r > 
0.3A) is 

(18.2 ± 5.0) (1 + z) 7 , 7 = 0.58 ± 0.50. 

The observed (^p) for Lyman-limit systems is 

~ 0.5 ±0.3 atz = 0.5 



~ 0.7 ± 0.2 with < z LLS >= 0.7 

(Stengler-Larrea et al. 1995; see also Mo & Miralda-Escude 
1996 for model predictions). The observed (^p) for damped 
Lya systems has been given above. The results of LBT W 
& CLWB show that at least 30 per cent of the absorbers 
with W r > 0.3 A are related with luminous galaxies. With 
comparison to these results, we summarize: (1) model B can 
be ruled out because the predicted (^p) of strong Lya lines 
are small and no damped systems can be predicted in such 
models. (2) model C can be ruled out for lack of damped 
systems and for insufficient (4^-) of Lyman-limits systems. 
(3) models D, F2 and F7 can be excluded because they give 
insufficient (^p) of strong Lya lines. (4) models B7 and F6 
predict too many Lyman-limit systems and should be ruled 
out also. Thus, we come to conclusion that models Fl, F3, 
F4, F5 are reasonable models, because they predict plausible 
number densities for strong Lya lines (W r > 0.3A), Lyman- 
limit systems, and DLAs. 



2. 6. 1 Velocity spread of absorption systems 

The line number densities above are predicted by mea- 
suring the overall REW of all the subcomponents of a 
galaxy /absorber. There are some cases where more than one 
subcomponent in a galaxy is crossed by the LOS. Some- 
times it happens that the velocity spread of a line system 
is very large. When the velocity spread of a line is less 
than 100 — 200 kms -1 , the HST key project would possibly 
have detected just one line, but otherwise this would have 
been counted as 2 lines. To study this effect, we plot the 
cumulative distribution of velocity spread of the simulated 
galaxy /absorber systems for model F3 in Fig.jj]. As we can 
see, about 78 (87) per cent of the systems have a velocity 
spread less than 100(200) kms^ 1 . So, if we treat those lines 
with velocity spread larger than 200 km s -1 as 2 lines, the 
predicted line number density will increase by about 15 per 
cent. 



o 1 

cd 0.9 

^> 
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Figure 1. The cumulative distribution of velocity spreads of the 
simulated galaxy /absorber systems for model F3. 



2.6.2 Photoionization flux contributed by galaxies 

Hot stars (O, B stars) in a galaxy can contribute some ion- 
izing photons (at A < 912A, namely Lyman continuum, 
hereafter Lyc), which may reduce the neutral fraction in 
the galactic halo. Here we make a simple estimation about 
the effect of extra photoionization by star formation in the 
galaxy itself. 

Stellar synthesis models suggest that the number of ion- 
izing photons emitted from a galaxy in a unit time (Q(H )) 
is related to the star formation rate (SFR) by 



SFR= 1.08 



Q(H°) 

10 53 s -l 



M 



©yr 



(see Kennicutt 1998 and references therein). Of the Lyc pho- 
tons, only ~ 14 per cent can escape the OB association and 
enter the diffuse ionized medium ("Reynolds layer") above 
the galaxy disk (Dove & Shull 1994b). About 65 per cent of 
the escaping photons are not absorbed in the HII layer and 
either escape from the galaxy or are absorbed by additional 
gas at high latitude (Dove & Shull 1994b). Thus approxi- 
mately 9 per cent of the Lyc photons can reach the outer 
halo. However, Leitherer et al. (1995) observed a sample of 
4 starburst galaxies and concluded that less than 3 per cent 
of the ionizing photons can escape. As an approximation, 
we use an escaping fraction f esc — 0.05. The rate of the es- 
caping Lyc photons is then fescQ(H°). The ionizing photon 
flux at galactic distance R is 

<b Lyc ~ 10 8 / esc Q5i#kpc cm ~ 2s ~\ 

where R kpc = fl/kpc, Q51 = Q(H°)/(W 51 s _1 ) and the 
energy flux P] is 



Jv — - — f&Lyc 
47T 



Lyc \ — 2 —1 — Itt —1 
1 ' 1 ' xlz , 



, 10' cm 2 s 1 / 
where h is Planck's constant. In dimensionless form we have 



J-21 — 10/escQsi^kpC 

For a normal spiral galaxy like the Milky Way Galaxy, 
with a SFR of 1 ~ 2M yr _1 , Q B i is about 100, and 
J_2i = 50i?^p C . At a distance of 10 kpc, J-21 = 0.5, which 
is larger than what we have used (J-21 = 0.05 at z = 0) by 
an order of magnitude. At a distance of 30 kpc, J_2i = 0.05, 
which is comparable to what we have used. At a distance of 



t If we take the shape of J v into account, the result should be 
multiplied by some factor. However because all the estimations 
here are quite uncertain, we omit this factor. 
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Table 2. Parameters and results of models 







N 


fl 


f2 


f3 


V dz Mean 


rdNs 

\ dz Mean 




Model 


parameters 


(total) 


(S) 


(SO) 


(E) 


(LLS a ) 


(>0.3A) 


fraction'' 


Bl 


M 5 = 1 


1357 


46.4 


24.1 


32.2 


0.40 


4.6 


21.1 ± 8.3% 


B2 


standard 


1078 


45.9 


22.3 


31.8 


0.45 


3.7 


17.1 ± 6.7% 


B3 


M 5 = 10 


925 


46.4 


21.6 


32.0 


0.35 


2.9 


13.4 ± 5.2% 


Bl 


Z = O.IZq 


844 


51.1 


19.3 


29.6 


0.42 


2.7 


12.4 ± 4.9% 


B5 


fg = 0.03 


807 


48.7 


21.2 


30.1 


0.20 


2.2 


10.1 ± 4.0% 


B6 


Vmi = 0.5 


1339 


47.2 


20.7 


32.1 


0.61 


4.6 


21.2 ± 8.3% 


B7 


J° 21 = 0.006 


1259 


51.3 


20.3 


28.4 


1.22 


4.9 


22.6 ± 8.9% 


B8 


SCDM 


472 


51.1 


16.5 


32.4 


0.16 


1.5 


6.9 ± 2.7% 


CI 


VDF1 


1984 


33.9 


23.3 


42.7 


0.07 


4.9 


22.6 ± 8.9% 


C2 


f g = 0.03,VDF1 


1897 


32.1 


26.3 


41.5 


0.06 


3.4 


15.7 ± 6.1% 


C3 


standard 


4454 


32.4 


24.3 


43.4 


0.15 


9.8 


45.2 ± 17.7% 


C4 


fg = 0.03 


4190 


32.0 


24.2 


43.8 


0.10 


6.1 


28.1 ± 11.0% 


Dl 


standard 


1041 


47.6 


20.4 


32.1 


0.48 


3.3 


15.2 ± 6.0% 


D2 


fg = 0.03 


741 


50.2 


20.2 


29.6 


0.36 


2.0 


9.2 ± 3.6% 


D3 


tUnf = 0.5 


1370 


44.0 


21.3 


34.7 


0.69 


4.9 


22.6 ± 8.9% 


Dl 


SCDM 


571 


50.4 


20.7 


28.9 


0.28 


2.1 


9.7 ± 3.8% 


Fl 


VDF1 


2591 


37.9 


23.1 


39.0 


0.70 


7.3 


33.6 ± 13.2% 


F2 


f g = 0.03, VDF1 


2335 


38.2 


21.9 


39.9 


0.44 


4.9 


22.6 ± 8.9% 


F3 


standard 


4885 


34.7 


23.5 


41.8 


0.69 


11.9 


54.8 ±21.5% 


Fl 


fg = 0.03 


4545 


34.3 


24.0 


41.7 


0.42 


7.7 


35.5 ± 13.9% 


F5 


Vint = 0.5 


5076 


33.8 


24.5 


41.6 


0.81 


12.9 


59.4 ± 23.2% 


F6 


J° 21 = 0.01 


5017 


35.4 


22.9 


41.7 


1.55 


20.1 


92.6 ± 36.3% 


F7 


SCDM 


3281 


33.3 


23.6 


43.2 


0.33 


6.1 


28.1 ± 11.0% 



a The Lyman-limit systems (LLS) are defined as absorbers with HI column density 10 17 cm 2 < Nm < 2 X 10 20 cm 2 . 
b The fraction of simulated lines related to galaxies to the number of lines per unit redshift at redshift z = 0.35, 
(dN/dz)(W r > 0.3A) ~ 21.7 ± 8.5 predicted by the HST QSO Absorption Line Key Project (Bahcall et al. 1996, see 
also CLWB). 



100 kpc J_2i = 0.005. For most of the halo clouds in our 
models, the typical distances to the galaxy center are larger 
than 30 kpc, so the additional photoionization flux by the 
galaxy itself may be negligible. This is also true for the satel- 
lites, because they are located at even larger distances. The 
typical absorption radius is about 100 kpc. If we make the 
extreme assumption that all clouds within 30 kpc are fully 
ionized, the number of lines would be reduced by about 10 
per cent. For early type galaxies, the SFR is lower and so 
the local photoionization can be neglected. The situation 
may be different for starburst galaxies with SFR bigger than 
~ lOM© yr _1 , but the number density of starburst galaxies 
is small and furthermore, the rate of emission of ionizing 
photons drops quickly in a short period of time (see Nulsen, 
Barcons, & Fabian 1998 for discussion). To estimate the ef- 
fect we have made simulations of the model F3 but with 
an ionizing flux increased by a factor of 2. The value of 
for strong lines is reduced to ~ 8.8 from 11.9. We therefore 
expect that our results will not be affected significantly by 
the local ionizing sources. 



2.7 Properties of absorbing galaxies 

Our simulations provide additional information about the 
galaxy /absorber systems. For example, Fig.^ presents the 
distribution histograms of the projected distances, equiva- 
lent line widths, luminosities, absolute magnitudes, circular 
velocities and redshifts of absorbing galaxies for the stan- 
dard model F3. The results are summarized as follows: 



(i) The REW distribution has a peak near 0.3A, the ma- 
jority of REWs are between O.lA and I.OA. 

(ii) About 70 per cent of projected distances are between 
20/i -1 kpc and 200ft" 1 kpc, only 10 per cent are at pro- 
jected distances less then 20ft" 1 kpc and 20 per cent are at 
projected distances larger then 200/i _1 kpc. 

(iii) About 70 per cent of the luminosities are between 
0.1-Lb* and l.OLs*, which implies absorbing galaxies are 
relatively luminous galaxies in the model. 

(iv) At least 80 per cent of the circular velocities are be- 
tween 100 km s" 1 to 300 km s _1 , the fraction of absorbers 
with circular velocities less than 100 km s" 1 is only about 
15 per cent. 

(v) The number density for LLS is about 0.69 which is 
in good agreement with observations (Stengler-Larrea et al. 
1995). 

(vi) The distribution with redshift is fiat and the average 
n(z) for absorbers with W r > 0.3A is about 11.9 which can 
account for about 55 per cent of the sources observed by the 
HST. 

Note that in the lowest right panel of Figj2|, the pre- 
dicted number density is almost independent of redshift. 
This result allows us to average the number density over 
the whole redshift interval in Table ^. In observation, the 
number density of strong Lyo? line evolves slowly with red- 
shift at low redshift (z < 1). It is valuable to investigate the 
correlations of REW versus projected distance, galaxy lu- 
minosity, galaxy /absorber redshift and study the fractions 
of absorbers produced by different morphological types of 
galaxies. 
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Figure 2. The distributions of equivalent line widths, projected distance, absolute magnitudes, circular velocity and rcdshifts, for model 
F3. In the lowest-right panel, the dotted line represents the results of the HST QSO Absorption Line Key Project (Bahcall et al. 1996), 
which is ^ = (18.2 ± 5.0)(1 + z) - 58 . The error bar is due to the uncertainty of the HST result. The number of lines per unit redshift 
for Lyman-limit systems is also shown as a dashed line. 



2.7.1 Projected distance 

We investigate the anti-correlation between REW, W r and 
projected distance of LOS to galaxy centre, p. A power-law 
relation is adopted, 

logW r = -alogp + C, (44) 

where a is the slope and C is constant. The results are sum- 
marized in Table |^. In the table, we list results for typical 
models B2, CI, C3, Dl, Fl, F3 (the standard model), and 
F5. In model Fs, we still get an apparent anti-correlation 
and a slope of ~ -0.4. Note that in all models the majority 
of absorption lines by halo clouds and/or satellites definitely 



decides the character of the anti-correlation (for example, see 
Fig.^) . The difference in the results found in r s and r v is rea- 
sonable because these three types of absorption components 
which do not necessarily have the same pattern of depen- 
dence. This difference has been noticed before (Le Brun et 
al. 1996). 

We plot the distribution of W r and p for some models. 
Panels in Fig.^ are for the results of model B2, CI, Dl, F3 
respectively. 

Small galaxies may have almost no cold gas because 
gas could be heated by supernova explosions or all the gas 
may have cooled to form stars, and there is no more gas 
to accrete. Thus these galaxies cannot contribute to absorp- 
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Figure 3. Example Model B2, CI, Dl, F3 (ACDM): The solid lines in both panels are for the linear fit logVK r = C — a log p. In each 
panel, about 95.4 per cent of the data points lie between the two dashed lines drawn with 2a c vertical shifts. The squares and circles are 
for spiral galaxies and clliptical/SO galaxies respectively. Small open symbols represent galaxies of luminosity Lb < O.lLs*, small filled 
symbols represent galaxies of luminosity O.ILb* < Lg < Lb*, large filled symbols represent galaxies of luminosity Lb > Lb*- 



tion. In practice, if we exclude those with V c h < 100 km s~ 
(as we can see in model F3, there are about 80 per cent of 
the absorbers with circular velocity of 100 km s _1 or more), 
there will be a stronger anti-correlation and a steeper slope 
of the linear fit at a highly significant level for almost all 
models. For model Fs, typically the Pearson co-efficient, r p 
is ~ -0.5 and the slope is ~ -0.5 at very high significance. 
This can be seen clearly from the results for 'sample b' in 
Table |. 

In summary, our predicted anti-correlation is significant 
but the dependence of W r on p is not as strong as that of 
CLWB. However, as pointed out by Tripp et al. (1998), se- 



lection effects could artificially tighten this anti-correlation. 
We will discuss this problem in the next section. 

2. 7. 2 Galaxy luminosity 

As suggested by CLWB, there is a power-law relationship 
between W and p and Lb, 

logW r = -alogp + /31ogL s + C. (45) 

We apply the analysis to model Fl, F3 and F5. The results 
are listed in Table ^. For example, for model Fl, the anal- 
ysis for sample b yields, C = 1.21, a = 0.79, /3 = 0.19. We 
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Table 3. The dependence of REW on projected distance 



Model 


c 




a 




r s (SL C ) 


r p 


(SL) 


c c 


B2 sample a a 


.29 ± 


.05 


.40 ± 


.03 


-.36 (11.8) 


-.42 


(13.8) 


0.26 


B2 sample b 6 


.66 ± 


.05 


.60 ± 


.03 


-.48 (15.4) 


-.57 


(19.8) 


0.23 


CI sample a 


.24 ± 


.04 


.37 ± 


.02 


-.33 (14.2) 


-.37 


(16.4) 


0.26 


C3 sample a 


.06 ± 


.03 


.27 ± 


.01 


- 28 (18 0) 


-.34 


(22.5) 


0.24 


Dl sample a 


.56 ± 


.06 


.56 ± 


.03 


-.35 (11.6) 


-.47 


(16.3) 


0.31 


Dl sample b 


.97 ± 


.06 


.77 ± 


.01 


-.45 (14.8) 


-.59 


(21.1) 


0.28 


Fl sample a 


.54 ± 


.04 


.51 ± 


.02 


-.42 (20.7) 


-.46 


(23.4) 


0.27 


Fl sample b 


.95 ± 


.05 


.70 ± 


.02 


-.49 (22.1) 


-.57 


(29.1) 


0.26 


F3 sample a 


.42 ± 


.03 


.42 ± 


.01 


-.38 (26.5) 


-.48 


(35.9) 


0.26 


F3 sample b 


.65 ± 


.03 


.52 ± 


.01 


-.43 (28.9) 


-.55 


(40.0) 


0.25 


F5 sample a 


.43 ± 


.02 


.43 ± 


.01 


-.44 (32.2) 


-.50 


(38.6) 


0.25 


F5 sample b 


.69 ± 


.03 


.53 ± 


.01 


-.50 (36.5) 


-.58 


(45.0) 


0.24 



a In sample a, we include those galaxy-absorber pairs with W r > 0.1 A. 

b In sample b, we exclude those galaxy-absorber pairs with circular velocity of central galaxies 
V cir < 100 km s" 1 and W r < 0.1A. 

c The statistical significance level, which is equal to r^J where r is co-efficient and TV 

is number of points. 



can determine the absorption radius of a galaxy r oc L B 
(t = P/a) with t ~ 0.24 - 0.32. The value of t is compa- 
rable but smaller than 0.37 which was derived by CLWB. 
However it is similar to that derived from Mgll obsorbers 
(Bergeron & Boisse 1991; Bergeron et al. 1992; Le Brun et 
al. 1993; Steidel 1995). Again, selection effects could lead to 
misleading conclusions (see discussion in 



1.3) 



2.7.3 A bsorber redshift 

We also analyse the dependence of the line equivalent width 
on absorber redshift assuming 



log W r = — alogp — 7log(l + z) + C, 
and 



log W r 



-a log p + /3 log L B — 7log(l + z) 



C. 



(46) 



(47) 



We apply the analysis to model Fl, F3 and F5. The results 
are listed in Table g and Table ^[ In summary, the relation- 
ship between REW and projected distance together with 
absorber redshift is marginally superior (with larger \r p \ or 
|r s |) to the relationship between REW and the projected 
distance but marginally inferior (with smaller \r p \ or \r 3 \) 
to the relationship between REW and projected distance 
accounting for Lb, and the anti-correlations between REW 
and projected distance accounting for Lb together with z 
is superior to the relationship between REW and projected 
distance accounting for Lb- For the analysis of eq.(k^), we 
have a ~ 0.48 - 0.80, j3 ~ 0.14 - 0.21, 7 ~ 0.49 - 0.61. This 
means the absorption radius of a galaxy r oc L B (1 + z)~ u 
(t = 0/a, u = 7/a) with t ~ 0.24-0.32 and it ~ 0.76-1.02. 
Our result of dependence on absorber redshift is different 
with the result of CLWB. CLWB concluded that REW do 
not depend on absorber redshift. However as we will discuss 
below, selection effects should be considered in the imaging 
and spectroscopic survey and the total redshift interval of 
LOS in observation may be not large enough to determine 
the relations. 



2. 7.4 Covering factor 

From eq. (^), and the results in Table [| we can estimate the 
average absorption radius and covering factor from redshift 
of to 1. For example, for 'sample a' (see notation in Table 
^|for its meaning) of model Fl, absorbers with REW larger 
than 0.3A follow the relation 



- = (— Y 

r, \LbJ 



(1 + z) 



(48) 



r 

r, \Lib* 

where r, = 0.23ft" 1 Mpc and t = 0.32, u = 0.86. Thus we 
can calculate the total number by integration 

1 /-i , \2-2it 

Motel = / nr, 2 C(1 t/s Ci(z)dz 



H{z) 

x 4>*T(l + 2t- s,LB min /L B *), 



(49) 



where Ci(z) is the covering factor and F is the incomplete 
gamma function. Note that for these absorbers with W r > 
0.3 A within r», it is not necessary that the covering factor 
is always larger than one, that is, for a LOS with a large pro- 
jected distance it is not always possible to find an absorber 
with W r > 0.3 A. If we choose s = 1.1, L Btnln /L B * = 0.007, 
0* = 0.027ft 3 Mpc" 3 then we get AT total = 20.4 x F(z), where 
= Ci(z)dz. Comparing with the total 



no = So 



(i+jO" 



number listed in Table Q, A^otai should be 7.3, the average 
value of F(z) should be about 0.36. However the integration 



r 1 

Jo 



(1+ Z )°- 28 



— -dz is about 0.85. Thus the average covering 

Va +«o(i+z) 3 

factor within 230ft" 1 kpc should be ~ 0.42 if Ci(z) can be 
treated roughly as a constant. The effective gas absorption 
radius should be \/0.42 x r* ~ 150ft" 1 kpc. Therefore the 
covering factor within 250ft" 1 kpc is (150/250) 2 ~ 0.36. 

Our predicted covering factor is in good agreement with 
the LBTW paper. In that paper, almost every galaxy with 
p < 70ft" 1 kpc gives rise to absorption, about five of 10 
galaxies with 70ft" 1 kpc < p < 160ft~ kpc give rise to ab- 
sorption, and just one of 9 galaxies with p > 160ft" 1 kpc 
give rise to absorption. Thus the covering factor within 
250ft" 1 kpc is about 0.31. However our predicted effective 
absorption radius (~ 150ft" 1 kpc) is a bit smaller than the 
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Table 4. REW dependence on projected distance and luminosity 



Model 




c 


a 




& 




r s (SL) 


r p (SL) 




Fl sample 


a 


.89 ± .04 


.63 ± 


.02 


.20 ± 


.01 


-.53 (27.8) 


-.59 (33.1) 


.25 


Fl sample 


b 


1.21 ± .05 


.79 ± 


.02 


.19 ± 


.01 


-.58 (29.7) 


-.65 (36.4) 


.24 


F3 sample 


a 


.62 ± .03 


.48 ± 


.01 


.15 ± 


.01 


-.46 (31.0) 


-.57 (44.4) 


.24 


F3 sample 


b 


.81 ± .03 


.57 ± 


.01 


.14 ± 


.01 


-.49 (34.4) 


-.61 (46.5) 


.23 


F5 sample 


a 


.66 ± .02 


.50 ± 


.01 


.15 ± 


.01 


-.51 (39.2) 


-.58 (47.4) 


.24 


F5 sample 


b 


.87 ± .03 


.59 ± 


.01 


.14 ± 


.01 


-.55 (41.9) 


-.63 (51.6) 


.23 



Table 5. REW dependence on projected distance and redshift 



Model 




c 




a 




7 




r 3 (SL) 


r v {SL) 




Fl sample 


a 


.62 ± .04 


.51 ± 


.02 


.16 ± 


.07 


-.44 (22.2) 


-.48 (24.4) 


.27 


Fl sample 


b 


1.03 ± 


.05 


.71 ± 


.02 


.49 ± 


.07 


-.51 (25.3) 


-.58 (30.2) 


.26 


F3 sample 


a 


.49 ± 


.03 


.42 ± 


.01 


.45 ± 


.05 


-.41 (29.2) 


-.50 (37.5) 


.25 


F3 sample 


b 


.74 ± 


.03 


.52 ± 


.01 


.48 ± 


.05 


-.46 (31.4) 


-.57 (42.1) 


.24 


F5 sample 


a 


.52 ± 


.02 


.43 ± 


.01 


.50 ± 


.04 


-.47 (35.1) 


-.52 (40.7) 


.25 


F5 sample 


b 


.78 ± 


.03 


.54 ± 


.01 


.50 ± 


.05 


-.53 (39.2) 


-.60 (47.0) 


.24 



174/i _1 kpc derived by CLWB. The reason of this difference 
could be the large covering factor used in CLWB. Indepen- 
dently, Bowen et al. (1996) derive a covering factor of 0.50 
within p < 160ft -1 kpc (three of six galaxies give rise to 
absorption) . 



out comparisons with results of imaging and spectroscopic 
surveys, it is absolutely necessary to construct absorber cat- 
alogues with observational selection criteria and investigate 
the possibility of mis-identification of absorbers (i.e., opti- 
cally unseen absorbers mis-matched with a bright neighbour 
galaxy). 



2.7.5 Galaxy morphological type 

CLWB conclude that galaxies which produce Lya absorp- 
tion systems span a wide range of morphological types from 
elliptical or SO galaxies through late-type spiral galaxies. 
Consistently the models show every type of galaxy can pro- 
duce Lya absorption systems. The predicted fractions of ab- 
sorbing galaxies in spiral, SO, elliptical galaxies are defined 
as fl, fl, f3 respectively and listed in Table |^. As we can see, 
the fractions vary from model to model. In model B and D, 
about half of the absorbers arise in spiral galaxies. In model 
C, this fraction is only about one third because E/S0 galax- 
ies possess more satellites than spirals. In model F, about 
35 per cent of the absorbers are produced by spirals. It is, 
however, not clear, if there are many clouds inside haloes 
of ellipticals/SO galaxies and how many absorbing systems 
can arise from tidal tails related with spirals. If we exclude 
absorption by halo clouds of ellipticals and SO galaxies, for 
instance, comparing models B2 and Fl, the fraction by spi- 
rals is estimated to be (7.3 x 38%)/(7.3 - 3.7 x 54%) ~ 52 
per cent. A crude estimate using Fig. 2 of CLWB shows that 
this number is about 70 per cent, however the information 
on galaxy morphology in observations may be inadequate to 
tackle this problem. We suggest that these numbers should 
be investigated further to see whether the gas absorption 
sections of different types of galaxies are the same or not. 



3 SELECTION EFFECTS: MOCK IMAGING 
AND SPECTROSCOPIC SURVEYS 

In section 2 we have studied (^p) an d the overall proper- 
ties of absorbers limited only by the lowest line width of 
O.lA, but without considering selection criteria. In carrying 



3.1 Selection criteria of absorber-galaxy pairs 

Our galaxy-absorber pair selection criteria are similar to 
those of CLWB. We only select absorbers with a rest frame 
line width 

W r > O.lA, (50) 

and with the B band luminosity satisfying 

m B < 24.3. (51) 

Similar to LBTW, we only select those galaxies within an- 
gular distances to the QSOs satisfying 

6 < 1'3. (52) 

A small velocity separation between absorber and galaxy 
centre, v = |c% a i — cz a bs| < 500 km s" 1 , is required to relate 
an absorber with a luminous galaxy (Lanzetta et al. 1997; 
CLWB; Bowen et al. 1996). This small velocity separation 
excludes almost all random galaxy-absorber pairs. 

3.2 Mis-identification 

LBTW argued that it is unlikely that the absorbing gas 
arises in dwarf companions of the luminous galaxies because 
no such dwarf galaxy was found in the LOS toward QSOs at 
redshifts z < 0.2 down to a luminosity of « 0.05L, [§]. But 

t See LBTW's Fig. 2. for the variation of projected distance 
threshold with redshift for an angular distance threshold 9 = 1'3. 
§ In our models, a spiral galaxy with Lb ~ 0.05Lb„ may have 
V v i r ss 55kms — 1 and mg m 23.2 at z ~ 0.2 in a ACDM cosmol- 
ogy, for example. 
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Table 6. REW dependence on projected distance, luminosity and redshift 



15 



Model 




c 




a 








7 


r s (SL) 


r p (SL) 


CT C 


Fl sample 


a 


1.01 ± 


.04 


.65 ± 


.02 


.21 ± 


.01 


.56 ±0.06 


-.55 (29.6) 


-.61 (34.7) 


.25 


Fl sample 


b 


1.34 ± 


.05 


.80 ± 


.02 


.20 ± 


.01 


.61 ±0.07 


-.60 (31.7) 


-.67 (38.3) 


.23 


F3 sample 


a 


.71 ± 


03 


.48 ± 


.01 


.15 ± 


.01 


.49 ±0.04 


-.49 (36.3) 


-.58 (46.3) 


.24 


F3 sample 


b 


.90 ± 


03 


.57 ± 


.01 


.14 ± 


.01 


.51 ±0.05 


-.52 (37.0) 


-.62 (48.5) 


.23 


F5 sample 


a 


.75 ± 


03 


.50 ± 


.01 


.15 ± 


.01 


.51 ±0.04 


-.53 (41.9) 


-.60 (49.6) 


.23 


F5 sample 


b 


.96 ± 


03 


.59 ± 


.01 


.14 ± 


.01 


.51 ±0.04 


-.58 (44.6) 


-.65 (53.8) 


.22 



van Gorkom et al. (1996) and Hoffman et al. (1998) have lo- 
cated faint dwarf galaxies at the reds hifts of a few low-z Lya 



lines. In addition, as we can see in §2.1, the lower limit in 
the luminosity function can be as low as ~ 0.007L*. Satel- 
lites are even fainter with typical circular velocity around 
30kms"\ Of course, at very low redshift, these faint satel- 
lites can be identified in imaging surveys with a large field 
of view. However, an angular threshold of 1'3 at z = 0.1 
means in general a distance of 100ft -1 kpc so that many of 
these satellites are excluded from pair samples because they 
typically have large distances to the central galaxy. 

As suggested by Tripp, Lu, & Savage (1998), the Lya 
lines could be due to undetected faint dwarf galaxies that 
are clustered with the observed luminous galaxies. This kind 
of Lya lines arise at small projected distances, but the host 
galaxies are too faint to be identified in galaxy imaging sur- 
veys, especially at modest to high redshifts. These un-seen 
absorbers (optically uncatalogued galaxies) can cause erro- 
neous identifications, since one could make a mistake to re- 
late the corresponding line with a nearby luminous galaxy 
at a larger projected distance. We simulate nearby luminous 
galaxies around a central galaxy by the two-point correlation 
function of normal galaxies, which is 



£(r) = 



5/jT 1 Mpc 



(53) 



where r is the separation of two galaxies and j ~ 1.8. Thus 
the galaxy number in a volume with radius of ft -1 Mpc is 



N(r, z) = n c (l + zf \ 4nr 2 [l + £(r, z)]dr, 



(54) 



where n c is the co-moving galaxy density described in § [2.1^ 
and z is the redshift of the central galaxy. The result of 
the Canada-France Redshift Survey shows a strong redshift 
evelution of the galaxy two-points correlation function which 
is 



£(r,z)=t(r)(l + z) 



(55) 



where e~0-2atz< 1.3 (Le Fevre et al. 1996). Because 
e is uncertain observationally, we use e = 1.5 for simplic- 
ity (Shepherd et al. 1997). Our results are not sensitive to 
the choice, because the redshift range covered is small. The 
luminosity distribution of these galaxies is consistent with 
the luminosity function and they are distributed around the 
central galaxy following f(r, z). The apparent magnitude 
of these galaxies is calculated from eq. (||). The distribu- 
tion of the pairwise velocity is a Gaussian with dispersion 
a = 400 km s" 1 (Efstathiou 1996; Mo, Jing & Borner 1993; 
Jing, Mo & Borner 1998). 

We define a galaxy- absorber pair intimately linking the 
absorber with a bright galaxy (an absorbing galaxy with 



apparent magnitude brighter than the luminosity limit, i.e. 
tub < 24.3) as a 'physically associated pair' or a 'physical 
pair' for simplicity. On the contrary a mis-matched galaxy- 
absorber pair is called 'spurious pair'. The method to find 
a mis-matched pair is as follows: When there is a faint ab- 
sorber (whose apparent magnitude is fainter than the lumi- 
nosity limit), its neighbours will be simulated to see whether 
there is a nearby bright galaxy (brighter than the luminos- 
ity limit) with larger projected distance (however not larger 
than 400/i _1 kpc). For a positive result, this bright neigh- 
bour will be chosen to pair with the absorption line arising 
from the faint absorber and the new projected distance will 
be chosen. For a negative result in the search of a bright 
neighbour, we classify the galaxy-absorber pair as a 'missing 
pair'. We define those bright 'physical pairs' and 'spurious 
pairs' as 'bright pairs'. We also call a pair outside a certain 
angular distance threshold as a 'missing pair'. For instance, 
at very low redshift some 'bright pairs' have large angular- 
separations to a QSO LOS. 



3.3 Effect on the W r 



Lb relations 



When comparing the W r — p — Lb relations for simulated 
galaxy/absorber pairs to the observations, the selection ef- 
fects mentioned above may have impacts on the statistics. 
For example, there could be some 'spurious' galaxies at large 
impact parameters within the redshift window of the ab- 
sorbers. In some of the surveys carrried out, the sky area 
surveyed is so large that there is always one galaxy (by 
chance) within the 500 kms -1 of the redshift. This may 
cause mis-identification of the absorbing galaxy and add 
noise to the correlations. On the other hand, faint absorb- 
ing galaxies without bright neighbours within the 500 km s - 
window will not be listed in the catalogues, which may re- 
duce the noise and strengthen the correlations. In observa- 
tions, the correlations are the results of the balance of these 
two effects. In principle, both effects may lead to misleading 
conclusions about the average galaxy/absorber distance. 

We simulate 200 sight lines as in section 2 using 
methods described above. For model Fl, our result shows 
that, if all bright galaxy/absorber pairs are used , we get 
logW r = 0.86 - 0.59 logp + 0.201og(L s /L s ») and X. = 

(rSl) 34 ' where r * = 220.7/i" 1 kpc (Wr > 0.3A). In 
contrast, for all physical pairs, we have log W r = 0.89 — 
0.631ogp + 0.201og(L s /L s «) and £ = (_^l) ' 32 , where 

r» = 174.8ft" 1 kpc (W r > 0.3A). As we can see, a larger r* 
is derived if spurious pairs at large distances are used. 
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3.4 Results for mock observations of known QSO 
LOSs 

We simulate observations for 10 QSOs at the redshifts listed 
in the CLWB paper. The redshifts of these QSOs are 0.200, 
0.264, 0.329, 0.371, 0.513, 0.534, 0.574, 0.616, 0.719, 0.927. 
The total redshift interval is about 5. Because of the proxim- 
ity effect, galaxies within 3000 km s _1 of the quasar redshift 
are excluded. With 100 mock observations, we get the sta- 
tistical number of galaxy-absorber pairs and the statistical 
properties of galaxy-absorber pairs. 

The total number of 'physical pairs' with W r > 0.3A 
for the ten QSOs is 

36.0 ± 6.3, 57.8 ± 7.5, 64.2 ± 7.0 

for model Fl, F3, F5 respectively. After applying the se- 
lection criteria, the total number of 'bright pairs' with 
W r > 0.3A is 

21.0 ± 4.8, 26.1 ± 4.8, 29.9 ± 5.3 

for model Fl, F3, F5 respectively. The predicted 'bright pair' 
numbers are in good agreement with that of CLWB paper 
in which there are 26 galaxies giving rise to absorption. 

If the model prediction is correct, we argue that the 
galaxy imaging survey at the faint end could be incomplete. 
This can be seen from the lower-left panel of Fig.^, in which 
the dotted line shows our predicted distribution of appar- 
ent magnitudes of 'bright pairs' and the solid line represents 
the distribution for all 'physical pairs', while the dashed line 
represents the distribution of apparent magnitudes of pairs 
from the CLWB paper. As we can see, considerable num- 
bers of the predicted ms of pairs are fainter than 25, and in 
the range of 23 to 25, our predicted 'bright pairs' are more 
numerous than the observed ones. In the lowest-right panel 
of Fig.^, the predicted number of 'physical pairs' in redshift 
bins at z < 0.5 is larger than that of CLWB, and compara- 
ble with that of CLWB at z > 0.5, while for 'bright pairs', 
the number is comparable with that of CLWB at z < 0.5 
but is lower than that of CLWB at z > 0.5. However the ob- 
served number density at high redshift is uncertain because 
the redshift interval there is quite small. Of course, if the 
difference is real, there are some implications. For example, 
the luminosity function at z > 0.5 would be higher than 
at low redshift and galaxies could be brighter in the blue 
band because of intense star formation. In conclusion, from 
the number of pairs and distribution of pair redshifts, the 
model predictions are consistent with current observations. 

In Fig.^, we also plot the distributions of p, W r , 
Lb/ Lb*, Mb as well as V c i r and compare the results with 
those of CLWB except for circular velocity. We apply the 
Chi- Square test to see if our predicted distributions are sim- 
ilar to those of CLWB (the data are from their table 4). The 
null hypothesis that the data sets are similar has probabil- 
ities of 97.0 per cent, 72.4 per cent, 65.7 per cent, 97.9 per 
cent, 31.1 per cent for W r , p, Lb/Lb*, Mb, z respectively. 

We give statistical results of the anti-correlation of W r 
versus p for all mock runs. Figj^ contains the histograms of 
the Spearman rank-order coefficient r 3 , the Pearson coeffi- 
ciet r v , and the zero point C, the slope a of the linear fit 
as in eq. (^). In the figure, the solid lines are for 'bright 
pairs' satisfying the selection criteria, the dashed lines are 
for 'bright pairs' with Kir > 100 krns" 1 , and the dotted lines 



are for all 'physical pairs'. Statistically the distributions of 
the slope a shift from small to large values for 'physical 
pairs', 'bright pairs', 'bright pairs' with Kir > lOOkrns -1 
and so on. Similarly, there are also shifts of lines in other 
panels. These shifts mean that selection effects can statisti- 
cally strengthen the anti-correlation of W r and p. We have to 
point out that this strengthening only has a statistical mean- 
ing and does not necessarily occur for every specific run in 
our simulations, because in some cases selection effects may 
also weaken the anti-correlations. Available results from ob- 
servations (CLWB; Tripp et al. 1998) are also shown with la 
error bars. Obviously some simulation runs can give consis- 
tent results compared with the observations. Note that the 
results of mock runs have considerable scatter, which may 
imply, as will be discussed later in that in the models 
the same total redshift interval as in present observations is 
not adequate to predict the real correlation. 

As mentioned above, the outcome of each run could 
differ case by case. Some examples of mock observations are 
given in Fig.|, Figg, Fig.|. 

The results for the anti-correlation of REW versus pro- 
jected distance are shown in Fig§ (for run No. 20). All the 
real galaxy-absorber pairs ('physical pairs') are drawn in 
the left panel. For the 108 'physical pairs' (of which 54 pairs 
have W r > 0.3A), we get the Spearman rank-order coef- 
ficient, r s — —0.41 (with significance level 4.6ct) and the 
Pearson coefficient r p = —0.62 (with significance level 8. Oct) 
and best fit line log W r = (0.85 ± 0.16) - (0.60 ± 0.07) log p. 
The galaxy-absorber pairs after applying selection criteria 
are drawn in the right panel. For 49 'bright pairs' (of which 
30 pairs have W r > 0.3A), we get r s = —0.58 (with signifi- 
cance level 4.9ct), r v = —0.75 (with significance level 7.7ct) 
and logWr = (1.63 ± 0.26) - (0.95 ± 0.12) log p. 

In comparison, CLWB give 

logW r = (1.34 ± 0.22) - (0.93 ±0.13) log p 

and Tripp et al. (1998) give 

logW r = (1.32 ± 0.20) - (0.80 ±0.10) log p 

by adding more LOSs. The results of 'bright pairs' for this 
run are in good agreement with those of CLWB and Tripp 
et al. (1998). As we can see, selection effects do strengthen 
the ant i- correlation in this specific run. However for some 
runs, selection effects do not strengthen the anti-correlation 
at all, because they only have a statistical meaning in the 
simulations. 

We analyse the relation between REW and galaxy lu- 
minosity using eq. ([IB]) for the run. The results are shown 
in Fig]?} In the left panel, C, a, fi are 0.96 ± 0.16, 0.63 ± 
0.07, 0.13 ± 0.04 respectively. In the right panel, they are 
1.80 ±0.26 of 1.02 ±0.12, 0.18 ±0.09 respectively. Again, in 
comparison, CLWB give 

log Wr = (1.78 ± 0.20) - (1.02 ± 0.12) log p 

± (0.37 ± 0.10) \og{L B /L B *)- 

As we can see, the zero point C and a of our prediction for 
'bright pairs' in run No. 20 are in good agreement with those 
of CLTW. However our results for f3 in this run are less than 
that of CLWB (still within the 2ct standard deviation), but 
in agreement with the value of 0.1 — 0.2 suggested by Bowen 
et al. (1996). 
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Figure 4. The distributions of equivalent line widths, projected distances, absolute magnitudes, circular velocities, apparent magnitudes, 
and redshifts for model F3. We compare the results with CLWB's except for the circular velocity distribution. For 'bright pairs', the 
distributions are in good agreement with those of CLWB. However, the numbers of apparent magnitudes for 'bright pairs' (dotted lines 
in the lowest left panel) in bins between 22.0 and 25.0 are higher than those of CLWB. This may mean that optical surveys could be 
incomplete in this apparent magnitude interval. 
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Figure 5. Statistical analysis of log W r = — alogp + C (for model F3). The upper left, upper right two panels are for distributions 
of constant C and slopes a for 'brigh t pairs' respectively. The lower left and lower right panels are for distributions of the Spearman 
rank order co-efficient r s and the Pearson co-efficient r p . All the lines here are only for pairs with W r > 0.1 A. The dotted lines are 
for all 'physical pairs'. The thick lines are for 'bright pairs' with 9 < 1'3. The dashed lines are for 'bright pairs' with 8 < 1'3 and 
V cir > 100 kms" 1 . 



There is a substantial number of missing pairs in our 
mock observations. Fig. || gives the results for the run. As 
we can see, there are a number of bright 'physical pairs' 
with angular distances to the LOS outside the threshold 
of 1'3 and some faint 'physical pairs' without bright neigh- 
bours located within 400ft -1 kpc. These pairs are defined as 
missing pairs and could not be listed in optical catalogues 
of absorbers. Note that there are also a few spurious pairs 
with angular distance outside the threshold. 



4 DISCUSSION 

As we remarked in the last section, if the total redshift in- 
terval is small, the results of correlations can have large sta- 
tistical deviations. In order to investigate this, we randomly 
choose 40 groups of LOS with the same number of LOS in 
each group from the simulation. We calculate the correla- 
tion coeffecient, the slope and zero point of the linear fit for 
every group. At last the average values over all the groups 
as well as the standard deviations can be calculated. Then 
we increase the number of LOSs in each group. These values 
will change until the total redshift interval is large enough to 
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get stable values with small deviations. We make plots of the 
dependence, from which we can determine how large a red- 
shift interval is adequate for an observational sample. The 
results for model F3 for 'bright pairs' are presented in Fig.^| 
Clearly for 'bright pairs' a total redshift interval ~ 10 is nec- 
essary to get statistically accurate value of r p , the Spearman 
coefficient, and a, the slope of the linear fit in eq. (144) , with 
standard deviation less than 0.1. If we want to determine 



the zero point (i.e. with standard deviation less than 0.2), 
the total redshift interval should be about 20 for model F3. 
This implies that results of the anti-correlation of W r versus 
p by different surveys may differ from each other if the total 
redshift interval in the survey is not large enough. 

In the models, only one third of the absorbers reside 
inside galactic haloes and two thirds of them are satelltes 
around central halos. This picture may reconcile different 
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Figure 9. The statistical analysis of the anti-correlation between 
W r and p (model F3) for 'bright pairs' with mj < 24.3 and < 
1'3. If the total redshift interval is 10, r p , C, a are —0.50 ± 0.11, 
0.55±0.26, 0.46±0.11 respectively and -0.54±0.11, 0.75±0.25, 
0.53 ± 0.11 respectively for the sample with V c i r > 100 kms -1 . 



conclusions by various authors (Morris et al. 1993; LBTW; 
CLWB; Bowen, Blades, & Pettini 1996; Le Brun et al. 1996; 
Tripp et al. 1998; Impey at al. 1999). Firstly, the models 
predict a large absorption radius and a reasonable covering 
factor. Secondly, the predicted absorbers are still closely re- 
lated with galaxies. Furthermore, it is possible that there are 
still a number of satelltes around central haloes even at dis- 
tance > 400/i -1 kpc. Thus if a catalogue of galaxy- absorber 



pairs includes those pairs with large projected distances, the 
anti-correlation of W r versus p could be weaken. 

As listed in Table ^[ our prefered models can explain up 
to ~ 55 per cent (and even more if an absorption system with 
large velocity spread is counted as 2 lines as would be the 
case with the HST spectral resolution) of the HST observed 
counterpart for lines wider than 0.3A. Can our models pre- 
dict more absorbers? Actually the evolution of the galaxy 
luminosity functions with redshift can change our predicted 
absorber number density. Because the numbers of aborbers 
in our simulations are directly proportional to the galaxy 
number density, our predicted number of absorbers can be 
higher at higher redshift if the galaxy number density there 
is higher. For example, if we chose the AUTOFIB luminosity 
function at 0.35 < z < 0.75 whose <f)* is about 1.5 times the 
(f>* at 0.02 < z < 0.15, the predicted number density could 
be as high as 17.8, which may account for 80 per cent of 
the observed counterparts. Furthermore, the higher galaxy 
density at higher redshift may also increase the absorber 
density at higher redshift in Fig.^| and solve the possible 
discrepancy in high-redshift absorber number between our 
prediction and that of CLWB. 

Caution should be taken with our results, not only be- 
cause there could be some alternative Lya absorption aris- 
ing in other parts of the galaxies, but also because it is un- 
clear under present circumstances whether there are many 
satellites in the vicinity of big central galaxies and whether 
they possess gas (see Klypin et al. 1999 for further discu- 
sion; see also Bullock, Kravtsov, & Weinberg 2000; Charl- 
ton, Churchill, & Rigby 2000). After all we did not include 
all possible absorption components related with galaxies. For 
example, Morris & van den Bergh (1994) suggested that a 
significant fraction of weak lines could arise from pressure- 
confined tidal debris in the enviroment of small groups and 
clusters of galaxies (cf. Mo 1994) . Tidal debris can increase 
the total gas cross section so as to increase the absorber 
number density and the corresponding absorption line can 
have a large projected distance of > lOOh" 1 kpc. However 
the absorption line number arising from tidal tails depends 
on the unknown gas cross section and the generally unknown 
lifetime of the gas in tidal tails. Thus it is not possible to de- 
temine exactly what fraction of absorbers arises in tidal de- 
bris. Another possibility is that some low surface brightness 
galaxies could possess huge gaseous haloes or discs which 
can also give rise to absorption. In addition, galactic wind 
is also another possible source (Wang 1995). Of course, the 
absorption by the IGM still could play an important role. 

In observations, it is difficult to assign a galaxy to an ab- 
sorber counterpart, because an imaging survey of galaxies is 
never quite complete down to the faint end. Also absorbing 
galaxies may be outside of the angular extent of the sur- 
vey. More LOSs to QSOs with higher resolution of the UV 
spectroscopy and more complete imaging surveys are nec- 
essary to investigate the physical origin and environment of 
the absorbers. At very low redshifts, it is possible to identify 
satellites with V c ir > 30kms _1 in optical imaging surveys, 
and then we can examine whether these satellites can give 
rise to Lya absorption lines. To discriminate models, it is 
also important to get more physical information about the 
absorbing components, such as size, temperature, metallic- 
ity, ionizing parameter, rather than only informations about 
W r and p. Observations of multi-LOS (QSO pairs or lensed 
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images) can provide useful tools to get more insight into 
absorbers (see Rauch 1998 and references therein). Further- 
more, the possible detection of line emission from extended 
gas in galactic haloes may al so help to determine the prop- 
erties of the gase ous haloes (Cirkovic, Bland- Hwathorn, & 
Samurcvic, 1999). 



5 SUMMARY AND CONCLUSIONS 

In this paper, we present results of Monte-Carlo simulations 
of Lya absorption line systems at redshift z < 1. To get 
constraints on the parameters, we simulate a set of models 
with different absorption components and various parame- 
ters. We compare the predicted absorption line densities for 
strong Lya lines, Lyman- limit systems as well as damped 
Lya systems with observations. From these comparisons, 
some models can be excluded. In summary: 

(i) Models with a single absorption component (galaxy 
disc, or cold clouds in a galactic halo, or satellites) or models 
with disk and clouds in a galactic halo, cannot explain the 
observed number densites for strong Lya lines, for Lyman- 
limit systems and for damped Lya systems at low redshift. 

(ii) Models with all three components (galaxy disc, cold 
clouds in a galactic halo, and satellites) can explain up to 60 
per cent of the observed number density for strong Lya lines 
at low redshift. These models can also predict reasonable 
number densities of Lyman-limit systems and damped Lya 
systems at low redshift. 

(iii) The fraction of the line number density for strong 
Lya lines due to satellites is ~ 40 per cent more than that 
due to clouds in galactic haloes (which is ~ 20 per cent) by 
about a factor of 2. The exponential galaxy discs can only 
account for a small amount of strong Lya lines. If indeed 
there are large numbers of satellites surrounding big central 
galaxies and they possess gas, these satellites may play an 
important role for strong Lya lines at low redshift. 

(iv) The predicted (ttt) for Lyman-limit systems due to 
cold clouds in galactic haloes is ~ 0.4, which can account for 
most of the observed Lyman-limit systems. The line number 
density of Lyman-limit systems due to satellites is only ~ 
0.1, which is four times smaller than that due to clouds in 
galactic haloes. 

The properties of the predicted absorbers, such as 
REW, projected distance, galaxy luminosity, circular veloc- 
ity and absorber redshift have been analysed. The predicted 
dependence of line width on projected distance is W r oc p~ a 
with a ~ 0.4 - 0.6, rather than a ~ 0.8 - 0.9 (cf. CLWB; 
Tripp et al. 1998). This predicted anti-correlation is weaker 
than the observed one because we include all faint absorbers 
(with apparent magnitude fainter than magnitude limit in 
optical surveys) which have small impact parameters and 
we include absorbers with impact parameters larger than 
200/i _1 kpc. Other correlations of REW versus luminosity 
and/or absorber redshift have also been investigated. In gen- 
eral, if we assume W r oc p~ a Lg(l + z)~~* , the analysis gives 
a ~ 0.5, j3 ~ 0.15, 7 ~ 0.5. This means the average absorp- 
tion radius of a galaxy r oc L l B (l + z)~ n with t ~ 0.3 and 
u ~ 1.0. The average covering factor within 250fo -1 kpc is 
estimated as ~ 0.36 which is in good agreement with pre- 
vious results (LBTW). The effective absorption radius is 



estimated to be 150/i _1 kpc, which is consistent with the 
observational result ~ 170/i _1 kpc derived by CLWB. 

To compare with results of imaging and spectroscopic 
surveys, it is necessary to study selection effects. Selection 
effects have impacts on the statistics of the galaxy/absorber 
properties. The present of 'spurious' galaxies at larger im- 
pact parameters within the redshift window of the absorber 
and the 'missing' of faint absorber at small impact param- 
eter may lead to misleading conclusions that the average 
galaxy/absorber separation is very large. Our simulations 
show that this is indeed the case. We construct mock obser- 
vations with the same known QSO LOS as CLWB applying 
selection criteria which are similar. By an adequate num- 
ber of mock runs, the total number of galaxy-absorber pairs 
can be predicted and the correlations mentioned above can 
be analysed. The predicted number of galaxy-absorber pairs 
with W r > 0.3A is ~ 26 ±5, in good agreement with CLWB 
(~ 26). The analysis of the anti-correlation between W r and 
p shows that selection effects can statistically strengthen 
the anti-correlation. Some results for mock runs can pro- 
duce anti-correlations consistent with CLWB. We also pre- 
dict some 'missing galaxy-absorber pairs' which are excluded 
by the selection criteria. 

We estimate the redshift interval adequate to predict 
accurate anti-correlations of W r versus p. To get results with 
a small scatter, it is found that in the standard model a total 
redshift interval of ~ 10 is required. This redshift interval is 
twice that of the LOSs in CLWB. This may imply that the 
total redshift interval in present surveys is not large enough 
to reveal the real anti-correlation. 
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APPENDIX A: SELF-SIMILAR SOLUTION FOR 
COOLING FLOW 



From eq. (18), we write, 



1 
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{x + u)x 2 dx 
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(x 2 u) 



p(x) 

(x + u)x 2 {l + x) 2 ' 
and get the general solution as 
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It is easy to solve this equation for the case that u does not 
depend on x. Then 



Pc 



8 



+ x)(b- l) a 



© 2000 RAS, MNRAS 000, 



22 W.P. Lin, G. Borner and H.J. Mo 



+ 



16 



24/n 



i+x 



Lanzctta K.M., Webb J.K. 



(b + k)(b-l) s (b + x) 2 (b~l) 2 (b-lY 
Here we let u = b. For x — > oo, p c (:r) — ■> 0, and have c = 0. 
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